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In this paper, we study quantum Sp(N) antiferromagnetic (AF) Heisenberg models 
in two dimensions (2D) by using the Schwinger-boson representation and the path- 
integral methods. An effective field theory, which is an extension of CP^^^ model 
in (2+l)D, is derived and its phase structure is studied by the 1/A^-expansion. We 
introduce a spatial anisotropy in the exchange couplings and show that the effective 
coupling constant in the CP^~^ model is an increasing function of the anisotropy. 
For the SU(N) AF Heisenberg model, which is a specific case of the Sp(N) model, 
we found that phase transition from the ordered "Neel state" to paramagnetic phase 
takes place as the anisotropy is increased. In the vicinity of the SU(N) symmetric 
point, this phase structure is retained. However as a parameter that controls explicit 
breaking of the SU(N) symmetry is increased, a new phase, which is similar to the 
spiral-spin phase with a nematic order in frustrated SU(2) spin systems, appears. 
It is shown that at that phase transition point, a local SU(2) gauge symmetry with 
composite SU(2) gauge field appears in the low-energy sector. It is another example 
of symmetry-enhancement phenomenon at low energies. We also introduce a lattice 
gauge-theoretical model, which is a counterpart of the effective field theory, and 
study its phase structure by means of the Monte-Carlo simulations. 



I. INTRODUCTION 



Study on quantum antiferromagnets has been one of the most active area in the condensed 
matter physics. In particular in the last decade, much attention has been paid to exotic phase 
and exotic phase transition for which Landau's classic paradigm cannot be applicable p]. It is 
expected that investigation of such exotic states is useful to understand anomalous properties 
of under doped high-T^ materials j^. Furthermore recent development in technologies of 
ultracold atoms and optical lattice trap elevates purely academic quantum spin models to 
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realistic ones, and these cold-atom systems are sometimes regarded as a final simulator 
for strongly-correlated electron systems. Quantum SU(N) antiferro magnets are one of these 
examples. Theoretically these models can be studied by using the Schwinger-boson methods 
and the 1/A^-expansion ^. In a seminal paper|^, it was shown that spin-(A^ — l)/2 (A^ is 
an even integer) cold atom systems in an optical lattice with one atom per quantum well 
can be regarded as quantum Sp(N) antiferromagnets. There are two parameters Ji, J2 in 
the Hamiltonian of the Sp(N) magnets and when Ji = J2, the symmetry is enhanced to 
SU(N)dSp(N) and a SU(N) quantum antiferromagnet is realized. 

In the present paper, we shall study Sp(N) quantum antiferromagnets by using the slave- 
boson (Schwinger boson) representation and the path-integral methods. We first derive an 
effective field theory for the Sp(N) Heisenberg model. This field theory is an extension 
of the CP^~^ model for the SU(N) antiferromagnetic (AF) Heisenberg model. Then we 
investigate its phase structure by using the 1/A^-expansion. We also study numerically its 
lattice-gauge-model counterpart by means of the Monte-Carlo (MC) simulations. 

This paper is organized as follows. In Sect. 2, we shall derive the effective field theory for 
the Sp(N) AF Heisenberg model by using CP^~^ representation of (pseudo-)spin degrees of 
freedom. In Sect. 3, we study the field-theory model by the 1/A^-expansion. We focus on 
quantum phase transition in the spatial two-dimensional (2D) system. We first show that 
in the SU(N) AF magnets with anisotropic exchange couplings on a square lattice, a phase 
transition from the AF Neel state to the paramagnetic state takes place as the anisotropy 
is increased. These two states persist in the Sp(N) system in the vicinity of the SU(N) 
symmetric point Ji = J2- As J2I J\ is decreased to some critical value, a phase transition to 
a new phase with a composite vector-field condensation takes place and a new kind of spin 
order appears. In Sect. 4, we study a lattice version of the obtained field theory and show the 
results of the numerical study for the Sp(4) case. Obtained phase diagram is qualitatively 
in agreement with that obtained by the 1/A^-expansion, but we also find some discrepancy 
between the 1/A^-expansion and the numerical study, e.g. order of the phase transition, etc. 
Section 5 is devoted for conclusion. In the Appendix, some details of the derivation of the 
effective field theory are given. 
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II. MODEL AND EFFECTIVE FIELD THEORY 

A. Sp(N) Heisenberg model 

We consider anisotropic Sp(N) Heisenberg model on a square lattice, i.e., exchange cou- 
plings in the x and y-directions are different. Quantum Hamiltonian of the system is given 
as follows in the most general form, 

^ = E {^1" Erf rf - $^r^r|}, (2.1) 

(ij) a,h a 

where i,j denote lattice sites and F"^ G sp(N) (the Lie algebra of Sp(N)) have ^^^"'"^^ 
components. On the other hand, G su(N)/sp(N), and there are (^+iK^~^) components. 
Hereafter in most of cases, we set the exchange coupling and J2'' as follows, 

4"+"^ = ^n,x, 4"^' = Jn,y, [u = 1,2), otherwise 0, 

where x (y) is the unit vector of the a;(t/)-direction. For the = 4 casejsj], the above 
generators are explicitly given as 

r« = (T"®/i^a = 1,2,3, F^ = l®/i^, F^ = l®/i^ (2.2) 

and 

r'' = i-[r^F''], (2.3) 

where a", fi"'{a = 1, 2, 3) are two sets of the 2x2 Pauli matrices. The Sp(4) case corresponds 
to the spin-3/2 system, and F° involves even powers of the SU(2) spin matrices, whereas 
-pafe involves odd powers of the spin matrices. Then the long-range order of F*^ indicates a 
spin-nematic order. 

It is useful to introduce the matrix J', which has the following properties and is a gener- 
alization of the time- reversal matrix ia'^ in the SU(2) spin case, 

J' = -J, = -1, JT^'^J = {T^y, JT'^J = -{ry. (2.4) 

For the Sp(4) case with Eqs.dMD, J = iay ^/j,"". 

We introduce the Schwinger boson operator ba {a = 1, ■ ■ ■ , N) and represent the "spin 
operators" in the Hamiltonian (12. ip in terms of them, f"'' = a^'^sbi^p, = blcJ'aB^hP- 
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By using the identities, 

y^^r^gF"^ = 25afj5p^ — IJa^J^a, (2.5) 

a.h 

r^^rjj^ = 2ba.a^p^ + IJa^J^f, — bafib^a^ (2.6) 

a 

we obtain 

7/ = Y.W^ - - 2(^i" + (2.7) 

where Qij = Ja^i,abj,p and K'ij = b^a ibaj- The operator Jai3bj,i3 is the conjugate spinor of 
6j,/3 and then Qij represents pairing of spins on the lattice sites i, j, whereas Kij corresponds 
to the Schwinger boson (spinon) hopping. Subsidiary condition 

N 

5^6j,,6,,Jphys) = |phys), (2.8) 

a=l 

must be imposed as the physical-state condition. We redefine the exchange couplings as 

J,,, ^ 2(Ji'^ + J^'^), ^ 2(jj'^ - J^'^), then 

From Eg. (12.91) . it is obvious that when J'-j = (Jjj = 0), the model have the global SU(N) 
symmetry. In the following section, we shall derive the effective field theory for the system 
(12. 9p by using the path-integral methods. 



B. Effective field theory 



In this subsection, we shall derive the low-energy effective field theory of the Hamiltonian 
(12. 9p . To this end, we use the coherent path-integral methods. For the Schwinger boson, we 
use the CP^~^ boson that satisfies Zi ■ Zi = ^i.a^i.a = 1 corresponding to the condition 
(12. 8p . Then the partition function is given as 

f/3 



Z = y T)zVz5{z ■ z — 1) exp 
M^) = - X] - 'H(z, z), 



dTA{T) 



(2.10) 



where 1-L{z, z) is obtained from Eg. (12. 9 p by replacing bi^a Zi^a and — )■ Zi^a- In order 





even site Q 
odd site A 



FIG. 1: CP^ variables on odd sites are paired with their NN z^. 

to obtain the effective field theory from Eq. fl2.10p . we integrate out the half of the CP^~^ 
variables, e.g. those at odd sites|6|]. To this end, we introduce a complete-orthogonal set 
of vectors in the CP^^"*^ space, {wf, i = 0,1, ■ ■ ■ , N — 1}. Arbitrary CP^""'^ variable z can 
be expanded as 2; = '^^PiWi, where p^'s are complex numbers satisfying = 1. As 

we are interested in the case with the both exchange couplings Jij, J[j > 0, dominant 
configurations in the path integral are given by Ze ^ z at even sites and Zo ~ J^z at odd 
sites, where 2; is a smoothly varying CP^^^ field. We parameterize odd-site Zo by referring 
to one of its nearest-neighbor (NN) even sites z^ as (see Fig{T]), 



(2.11) 



where Wq = Ze, Wi = J z^ and we do not have to specify the other Wj's. Then by the above 
remarks on the dominant configurations, we can assume pi(i 7^ 1) ^ 1 and 



pi 



f2.12) 



where ?7 is a U(l) variable that appears as a result of the local U(l) symmetry of the system. 



id 



From Eq.([2ll2]), 
By substituting Eg. (12. 131) into the action (I2.10p . we obtain 



.(i-iE,.,}.. 



A[r) 



(Ao + p-k + T-p-pMp) + 0(p^ 



(2.13) 



(2.14) 



odd— site 
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where p = {po,P2,P3, ■ ■ ■ iVn-iY ■ The matrix M is exphcitly given as follows, 

/ -4 - 4 j;^. • • • \ 



M 







-4 7 • 








V 



(2.15) 





■ ■ ■ -4 Jo,,- / 

and the "vectors" k and 1 have rather complicated form of composite of 2;, z, w and w, 
which are explicitly shown in the Appendix. In Eq. fl2.15p . Jo,j {J'o,j) denotes Jjj- (Jj'j) for 
i = o and the NN even j (see Fig 12]). As the kernel M has negative diagonal elements, the 
Gaussian integral over p can be safely done. 



Z = j VzVz5{z ■ 2 - 1) exp drA^T) 



VpVp exp 



rfrA(r) 



exp 



r r 

.Jo 



A.ir)= T(r)M-i(r)k(r)+ Ao{r) 



odd— site 



odd— site 



Ao = -Y^ Jo,j{ZjJ Wi){ZjJ Wi) - J'^j{wiZj){ZjWi] 



(2.16) 
(2.17) 
(2.18) 



where ^ • in Eq.f l2.18p denotes the summation over even site j around odd site o (see Figj2]). 



0^=3 



-o- 

J=2 



even site Q 
odd site A 



FIG. 2: Zo and its NN z,'s on even sites. 



We can consider a continuum limit of the effective model f l2.17p . as the remaining variables 
Zj's can be regarded as a smoothly varying field z{r) [tq = r, ri = x,r2 = y). Hereafter we 
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explicitly set the exchange couplings between adjacent spins as follows, 

Jx = Joi Jy = ^Joi J'x = J'y = ^J'o^ (2-19) 

where A is the parameter for the anisotropy between x and y directions. We take the limit 
/9 = |; — )■ cxD, introduce explicitly the lattice spacing a and rescale the coordinates as, 



r ^■ 



V2(l + A)Joar = cr 



a;^^/i±Aa; (2.20) 



where c is the "speed of light" of the system. Then we obtain action of the effective field 
theory for anisotropic Sp(N) AF magnets in two dimensions, 

S=Ygj dh[D^zD^z-^{zJD^z){zJD^z) + a{\z\^ -I)), (2.21) 

where Z)^ = — zd^^z, a is the Lagrange multiplier for the CP^^^ constraint, 7 is an 
increasing function of J'^/ Jq, 7 = J'q/ Jo + 0{{Jq/ JqY), and the coupling constant g is given 
by 

1 + A , , 

a. (2.22) 



From Eq.f l2.22]) . it is obvious that the effective coupling g has the minimum at the isotropic 
point A = 1. As we see in the following section, this means that the anisotropy tends to 
break the AF order of the ground state. 

As it is shown in the Appendix, there also exists a Berry phase in the action. 



= "^VTTx(^ " f ) / ^''-''^MDuzD.z). (2.23) 

The Berry phase 02.231) with the fractional coefficient depending on the anisotropy A and 
J'-coupling does not suppress effects of instanton in contrast to that with an integer co- 
efficient, i.e., Sb does not give any substantial effects on the phase structure and critical 
behavioral . For the case of the SU(2) AF magnets on 2D lattice with anisotropic couplings, 
this observation has been verified directly by the numerical study Q]. 

In the following sections, we shall study the field theory f l2.2ip by analytical and numerical 
methods. 
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III. PHASE STRUCTURE: ANALYTICAL STUDY 

A. expansion: Case of small 7 

The partition function of the effective field theory is given as, 

Z = JvzVzVaexp(^~^ J dh[D^zD^z-^{zJD^z)izJD^z) + - f )]) , (3.1) 

where we have introduced^the factor in front of the action to perform the 1/A^-expansion 
in the analytical study[ll]. At 7 = 0, the system (13. ip has the global SU(N) symmetry, 
z{r) — )• Vz{r), V G SU(N). We first consider the case of small 7, and put the following 
parameterization for z, z = zq+u+w, where zq = (no, 0, ■ ■ ■ , 0), u = {0,U2, ■ ■ ■ ,un) and v = 
(0, f2, ■ ■ ■ , f at). The fields u and v are real vectors. We have used the gauge-fixing condition 
vi = 0. After substituting the above parameterization to Eq. (l3.ip and keeping the quadratic 
terms, we perform the Gaussian integration over u and v as the leading order of 

N 
2^ 

^ VnoVaexp{-Ses{no,a)), (3.2) 
where 

TV /■ 

5'eff(no, a) = {N- l)Tr \og{-dl + a) + — j d'xa{nl - 1). (3.3) 

As the 7-term generates only higher order terms of u and v, it does not give any effect in 
the leading order of 

Gap equations are obtained as follows from Ses{no,<j) in Eq. (l3.3p . 



'Dn^'Da'DuDv exp 



dh{{d^uf + {d^vf + a{nl + + - 1)} 



6Sos{no,cT) N 

r = — ano = 0. (3.5) 

dno g 

We use the Pauli-Villars regularization with a cutoff A (~ vr/a) for the integral (13. 4p . and 
obtain the critical coupling gc by putting a = no = 0, 

- = —A. (3.6) 
9c 271 ^ ^ 



There are two phases, one for g > g, 

27r 

\g, g 



^ 27r / 1 1\ , ^ , ^ 

= ' ^0 = 0, ao = (a), (3.7) 
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and the other for g < gc, 



nl = l- —, (To = 0. 



(3.8) 



For g < gc, Sp(N) symmetry is spontaneously broken and 

J2{{z{r)r''^z{r)){z{r')r^^z{r'))) ^ 0, for |r - r'\ oo, 

a,b 

J2{iz{r)T^z{r)){z{r')T^z{r'))) ^ 0, for |r - r'\ oo. (3.9) 

a 

In later section, the above result will be verified by the numerical study of the lattice model 
for the effective field theory 9| . 



B. Case 7 1: Auxiliary fields 

In this subsection, we shall consider the case 7 «i 1. It is useful to introduce two kinds of 
auxiliary vector fields and to investigate the phase structure of the model. Inserting the 
following identities to the partition function (13. ip (where Zi, Z2 are irrelevant normalization 
constants and will be ignored hereafter), 

AT r ^ 



1 = Zi J VXf, exp 
1 = Z2 VLdDLoexp 



N r 

— J dh{\^ - izd^zf 



29 



£r\u^ - zJD^zl"^ 



we obtain, 



Z = j VzVzVaVX^VQ^Voj^ exp ^ j £r[z{-dl + iA^^^ + a)z 



z{-fu^Jd^)z - z{'^Uf,Jd^)z + A^ + 7|w 



= j VzVzVaVX^Vcu^Vcu^exp ^ j d'^rlzi-dl + iA^V^ + AJ 
+ 7|wmP + ~ z{-fUf,Jdf,)z - z{-fUJf,Jdf,)z - cr] ) , 



where zXf, d f,z = \f,{z ■ d^z - d^z ■ z). 



(3.10) 
(3.11) 



(3.12) 



1. Strong- coupling region 

First we shall study the model (I3.12p in the strong-coupling region g > gc, in which 
hq = 0, {(t) > 0. We put B = •yUfj^J'd^, B = •yu^J'dfj, for notational simphcity. Then the 
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partition function is expressed as 



+ a + 7|w^p + B + B)u + u{-\^ d ^ + i{B - B)}v 



+ v{X^ d ^ + i{B- B)}u + v{-d'^ + + a + 7| 



-{B + B)}v-a\ 



(3.13) 



where z = u + iv. As the action in Eg. (13. 131) is a quadratic form of z = u + iv, integration 
over z can be done. We define 



/9 = V(^.K)' + 7^K'^ 



2 

Ml ! 



(3.14) 



and the resuh after the integration over z is given as 

/./V ./V /" 

VaVX^Vco^Vco^ exp ( - y Tr log(a2 + (3^) + — j dx^a^ , (3.15) 



where 



log(a2 + /32) = log {(-9; + + 7|c^,r + + (Am'^m)' + i'K'^.I'}- 



From Eq. (l3.15p . we shall obtain an effective action of the vector fields A^ and o;^ in power 
of them. To this end, we use the identity 



log(a2 + /32) = 21og(-9^ + a) 



+ Tr 
+ Tr 



27 



r 



-dl + a) 



I ^ |2 



+ 0(Aj,u;;,Aju;J), 



(3.16) 



for positive constant a. We explicitly evaluate the momentum integrals in Eq. (l3.16p . 

2,0 1 ^ . 



Tr 



-A,+ 



(3.17) 
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Tr 



27 



7 



(27r) 



Tr 



7^ 



^ 12 



where 
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(27r) 



(27r)3 

(P + 2g)^(p + 2g), 

g2 + a 2((p + g)2 + a)(g2 + a) 



(3.18) 



(3.19) 



5, 



'^.^-^)n(p), 

n(p) = ^(p^ + 4a)E(p)-^, 



(3.20) 



E(p) 



1 . -1 /P 
——-tan \/— , 

47r p V 4(7 



111/ 



(27r)3g2 + (x 

rf^g (p+2q)^{p + 2q). 



(27r)3 2((p + g)2 + a)(g2 + a)- ^^'^^^ 
Then effective action of the vector fields A^^^ and o;^ is obtained as follows up to the quadratic 
order of them, 



si^{X„uj,) = Nj ^A,(p)n^.(p)A.(-p) 



+ A^7 



2 / d^P .-. 



(2vr)^ 



uj^{p)Uf,^{p)uj^{p) 
cP'p 



+ N^{1 - 7) / -^c.,(p)Q(^)a;.(p) 



AT 



d^p 
2gJ ("2^ 



(7. 



For positive a and p <C 1, 



S(p) 

n(p) 
n,..(p) 



1 



47rp V 4cr 

l(p2 + 4^)^^_ 
2^^ ^SttVo" 47r leTrVo" 



IGtt^/o" 

and therefore the quadratic term of cu^ becomes as follows at low momentum, 



(3.22) 

(3.23) 
(3.24) 
(3.25) 



N 



d^p _ 



IGttv^ 



^ni^^vip)- (3.26) 
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From Eq.f l3.26|) . it is obvious that behaves hke a massive vector field for 7 < 1, for 
7 = 1 it becomes massless and behaves hke a kind of gauge field, and finally for 7 > 1 its 
nonvanishing condensation is expected to occur. More detailed study on the case 7 = 1 will 
be given in later section, and it is shown there that a SU(2) gauge model reaUy appears. 

Let us study the case 7 > 1 somewhat in detail for the large- limit. Condensation of uj^ 
apparently breaks the rotational symmetry of the space (or ^-rotation of the square lattice) 



and also the U(l) gauge symmetry to Z2[l2|]. Here we assume (w^) = w 7^ without loss of 
generality. We also assume that u is real by the gauge symmetry of the system. Then the 
action of z{r) is given as 



29 



j d^r[-zdlz - -fuj{zjd^z) - -fuj{zjd^z) + a(|2p - 1)]. (3.27) 



The action 5*^ in Eq. fl3.27p can be diagonalized by introducing field ^(r) as 

z{r) = i={e-^--e(r) + e^^-^(z ^^(r))}, (3.28) 

Si = S-. = Yj d\{l{-dl + 7(1 - l)uj')i + a{U - 1)}. (3.29) 

From Eq.f i3.29p . it is obvious that the field ^(r) acquires its mass squared a' = (7+7(1 — 7)0;^. 

From 5*^ in Eq.f l3.29p . we can have a gap equation and determine the critical coupling gc. 
By integrating out ^(r), we obtain 

and 

S'Ja') f d^k 1 N , , 

""'^ - N j^.TT—J - W- = 0- (3-31) 



5a' J (27r)3p + (T' 2g 

Solution to Eq. fl3.3ip is obtained as 



V = 27r ( — - i ) (3.32) 



9c 9 
1 _ A 

9c 27r' 

The above critical value Qc is the same with that obtained for the case of small 7. 
In the large- A^ limit, spin correlations are obtained from Eq. (l3.28p . see Figj3l 



(3.33) 
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FIG. 3: (a) Correlation of T"^ (b) Correlation of 



2. Weak-coupling region 



We shall consider the weak coupling region g < Qc, in which the condensation of z{r) 
occurs. We assume a smooth second-order phase transition from the phase of u; = to that 
of w 7^ as 7 is increased, and estimate the critical coupling 7c[15|l- From Eq. fl3.28p . we 
divide z as z = zq + z, where Zq = (zi, 0, ■ ■ ■ ,0, zn), z = {0, Z2, ■ ■ ■ , z^-i, 0). Furthermore 
we put (cr) = 0. Substituting the above expression of z into the action, we obtain 

N 



S 



2^ 



zi-df, + Xf^ + ^\u^\' + a)z 

iX^z 'd' ^z - Q^zjd^z - u^zjd^z - a 
S + Sq, 



(3.34) 



where 



^i-dl + \l + l\uJ>.\^ + cr)z 



29 J 

+ iX^z d ^z- uj^zjd^z - u^zjd^z - a 



S, = - I d\ 



2g 



zo(-<9^ + A^ + 7|w^|^ + a)zo 



+ iX^ZQ d hZq — uj^zoJO^zq — uj^z^JO^Zq 



(3.35) 



(3.36) 
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Integration over Zr^ can be performed as in the previous case, 



2^52^2; exp(— 5) 



exp 



rfV (-9^ + A^ + 7|a;Mr + a) 



2<7 



For 5*0, we substitute the following expression for zq, 



z! = {zi,0, ■ ■ ■ ,0,zn), 



(3.37) 



(3.38) 



and obtain 



z'i-dl + Xl + ^\co,\' + a)z' 



+ iX^z' d ^z' - bj^z' Jd^z' - bj^z! Jd^z' 

+ 'S'int, 



(3.39) 



where 



5i 



int 



(A^ + 7|u;^|2 + a){iVo' + A^o(^' + ^')} 



+ i\Nod^{z! - z!) - -lu^NoJd^z' - ^u^NoJd^z' 



(3.40) 



We cannot integrate out z' exactly and therefore treat the cubic terms in S'int as a "perturba- 
tion" of the l/A^-expansion[ll]. Then the final expression of the effective action ^^^(A^, w^) 
is given as follows up the quadratic terms of A^, a and u^, 



p2 



A, 



IJ.U n 

7 



7 

1-7 

7 



I '^/i 



(3.41) 



where o;^ = + io;^. 



15 



It is obvious that for no 7^ case, only the transverse modes of (i.e., J^fiPf^'^f^ — 
0) can survive at low momentum. Furthermore its mass is renormalized by no 7^ as 
7[nQ + (1 — Then for 7 > 7c = 1 + 'Cl^); instead of z{r), condenses and 

nontrivial correlations of the "spin operators", T""^, appear as a result of (u^) 7^ 0. In this 
phase, 

{T''\r)) = (^(r)rt(r)) 

= nf cos(27wx) + nf sin(27wx) , (3.42) 

where nf = Re[{OJT''\0]^ nf = lm[{OJT^^{0]- On the other hand, (f'^(r)) = {0^^{0- 

IV. SU(2) GAUGE THEORY AT 7 = 1 

In the previous section, we found that the composite vector field o;^ behaves like a massless 
gauge field at 7 = 1, and for g < gc it acquires mass squared propotional to ng = ■ 2;) as a 
result of the Higgs mechanism. In this section, we shall explicitly show that the three real 
vector fields (A^, uj^, ujj form a SU(2) gauge field minimally coupled with at 7 = 1. This 
is another example of the symmetry-enhancement phenomenon, i.e., emergent symmetry at 
low energies. 

We start with the action 5^=1 in Eq.f l3.12l) . 

K-^l + (9 ^ + + \u}^\^ + a)z - z{fl^Jd^)z - z{u}^Jd^)z - a , 

(4.1) 

where we have put 7 = 1. Hereafter we explicitly consider the CP'^ case but generalization 
to an arbitrary is straightforward. We first redefine the CP'^ field Z{r) from the original 
z{r) as follows, Z{r) = {zi{r),Z2{r),Z4^{r),—Z3{r)y. It is easily verified that Z{r) is a CP^ 
field, Eli\Z^{r)\' = l. 

It is straightfward to verify the following equation, 

z'^^z = ZJ^s'^^Z, (4.2) 



S. 



7=1 



2^ 
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where 



/ 1 \ 
10 
0-10 
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Similarly 



oo^.{zjd^)z + uj^{zjd^)z = iu;^{ZJ:2 'd' ^Z) + lu^iZT.^ 'd' ^Z), 



(4.4) 



where 



J 



/ l\ 

0-10 
10 

\ -1 oy 











i 0\ 










i 




— i 








V 





—i 


0/ 



Si 



(4.5) 
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It is obvious that Sj(i = 1, 2, 3) satisfy the SU(2) algebra. Let us define SU(2) gauge field A^j, 
as = (A^, A^, A^) = (w^jO;^, A^), then the Lagrangian (14. ip can be rewritten as follows, 

N 



7=1 



2^? 



\{d^ + tJ:-A^.)Z\' + a\Z\' -a 



(4.6) 



Lattice gauge model corresponding to the above SU(2) gauge theory (14. 6 p is under study 
and result will be reported in a future publication. However phase structure of the system 
can be inferred by qualitative discussion. As in the usual CP^^^ model coupled with the 
U(l) gauge field, there exists a phase transition that separates ordered and disordered phases. 
However, as the SU(2) gauge field fiuctuates the hopping of the spinon Z{r) more strongly 
than the U(l) gauge field, the critical coupling gc{'y = 1) is expected to be smaller than 
gc{l = 0). From this consideration, we expect that the critical coupling gdl) is a decreasing 
function of 7 for 7 < 7c, though in the previous discussion for small 7 ^ 1 by the 1/A^- 
expansion we do not find any 7 dependence of Qc- As 7 exceeds 7c, the condensation of 
Lo^ tends to occur and fiuctuations of are suppressed. Moreover the original U(l) gauge 
symmetry reduces to Z2 gauge symmetry by the Anderson-Higgs mechanism and fiuctuations 
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of are also suppressed. Then gd'j) starts to increase at 7 = 7c. The above expectation 
will be confirmed by the numerical study of lattice-gauge model in the following section. 
Expected phase diagram is show in the (7 — 5')— plane in Fig HI 
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FIG. 4: Phase diagram in (7 — g) and (7/5 — (7)— planes. 



In the following section, we shall introduce a lattice model for the effective field theory 
with general value of 7, and study it by means of MC simulations. 



even site Q 
odd site A 



FIG. 5: Interaction terms in A/, in Eq. (j5.ip . 



V. NUMERICAL STUDY 



In this section we formulate the effective field theory (13. ip on a cubic lattice, and investi- 
gate its phase structure by means of the MC simulations. We exphcitly consider the Sp(4) 
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FIG. 6: (a)Enegry ii^ as a function of ci for C2 = C3 = 0. There exists a discontinuity at ci = 4.5. 
(b) "Specific heat" C as a function of ci for 02 = 03 = 0. There exists a sharp peak at ci = 4.5. 
System size L = 20. 

model. Action of the lattice model is given as follows, 

Al = Ci ^ ] Zr+ij,Urfj,Zr + C.C. 
+ C2 



(5.1) 



where r denotes the cubic lattice site, /x = (0, 1, 2) is the direction index and it also denotes 
the unit vector in the /^-direction. Field CP^ variables and f/r/^ is a U(l) gauge field 

defined on link {r,jj?), Ur^ ~ g^^MC*"). The parameters Ci oc 1/g, C3 oc 7/(7 and the C2-term is 
the lattice Maxwell term (the so-called Wilson term) corresponding to {d^Aj^ — dyA^Y in the 



continuum spacetime 



151]. We also added the C4-term on the diagonal lines in the 2D layers 



that represents the exchange couplings between spins on even sites of the original lattice. 
See FigJHl The partition function Zl is given as 

[Df/][D^D^]cpe^^ (5.2) 

where [DzDz\cp denotes the integration over CP'^ variables, and we use a specific parame- 
terization for them in the MC simulations. For the MC simulations, we used the standard 
Metropolis algorithm of local update. The typical statistics was 10^ MC steps per sample, 
and the averages and errors were estimated over ten samples. The typical acceptance ratio 
was about 50%. We also used multi-histogram methods to obtain reliable results near the 



phase transition point 



1^. 
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FIG. 7: Enegry distribution N[E] for the pure model. At ci = 4.498, N{E) has double-peak 
shape, whereas it has single peak at ci = 4.496 and 4.500. System size L = 20. 

We first consider the case of the pure CP^ model with C2 = C3 = C4 = 0, which corresponds 
to the anisotropic SU(4) AF magnet 10(]. We calculate the "internal energy" E = {Al)/L^ 



and the "specific heat" C = {{Al — E)"^) / to study phase structure, where is the lattice 
size and we impose the periodic boundary condition in moose of calculations. In FigJHl we 
show E and C as a function of ci. It is obvious that E has a discontinuity at ci ~ 4.5 and 
C has a very large peak at ci ^ 4.5, which indicates a first-order phase transition. In order 
to verify this observation, we calculated density of states N[E] that is defined as 



N[E] = j [DU][DzDz\cp 5{Al - E) e^K (5.3) 

Result in FigJT] shows that N[E] has a double-peak shape at ci = 4.498, whereas it has 
a single peak at the other couplings. This confirms the existence of the first-order phase 
transition in the CP^ model, though corresponding phase transitions in the CP^ and CP^ 
models are of second order. Study of the correlation functions of the "spin" operators given 



20 




FIG. 8: (a)C for C2 = 2.0, C3 = 0. System size is L = 16, 20, 30. (b)Scaling function <j){x) in 
finite-size scaling (j5.4p . Critical exponents and critical coupling are estimated as z/ = 0.8, a = 0.11 
and cic = 2.23. 

later on shows that phase transition from "Neel" to paramagnetic states takes place at 
ci = 4.498. 

QpAf-i jjiodel in the 3D spacetime was studied by the 1/A^-expansion and it was suggested 
that there existed a second-order phase transition from ordered to disordered phases as 
the coupling constant is increased. However the present investigation by means of the MC 
simulations shows that the order of the phase transition varies as a function of the parameter 
N. Similar phenomenon was recently observed some related model, e.g. multi-Higgs U(l) 
gauge model in 3D. We also studied finite but small C2 cases and found that the phase 
transition is still of first-order. However at intermediate value of C2, C exhibits the finite- 



size scaling (see FigJS]), i.e., the data o: 



follows with a scaling function 0(a:)[l6|, 



C(ci, L) for system size L and C2 = 2.0 can be fit as 



C(ci, L) = L<^/^</)(Li/^e), e = (ci - ci,)/cie, (5.4) 

where Cic is the critical coupling at infinite system size and estimated as Cic = 2.23. This 
fact means that the phase transition becomes of second-order as the value of C2 is increased. 

Let us turn on the ca-term and see how location of the phase transition varies. We studied 
the system by varying value of ci with fixed C3, and found clear signal of phase transitions, 
see, e.g. C for C2 = 2.0, C3 = 1.0 in FiglHl We also investigated phase structure of the 
system with ci fixed and C3 varied, and found that there is another phase transition line. 
See FigJi 
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FIG. 9: (a)C for C2 = 2.0, C3 = 1.0. System size is L = 14, 20, 30. Critical coupling is estimated 
at ci = 2.35. (b)C for ci = 3.0, C2 = 2.0. System size is L = 18. There are two phase transitions 
at C3 = 3.7 and 5.2. 




1.0 2.0 3.0 4.0 6.0 C'j 




FIG. 10: Phase diagram in (cs — — ) plane for C2 = 2.0. Solid lines are phase transition lines 
obtained by measurement of E and C. Dots denote phase transition points actually observed by 
measurement. "Spin" correlation functions exhibit different behavior in regions I ~ VI. Similar 
phase structure is obtained for C2 = 0. 

Obtained phase diagram in the (cs — ^) plane is shown in FigJTOl There are four phases, 
which are identified by the measurement of E and C. We also investigated the behavior of 
the correlation function of the "spin" operators, 

G^{r) = J2{^'+r'^r'), G,{r) = Y,{^X'^^^)- (5-5) 

a ab 

As we expect that a phase transition to a "spiral state" takes place as the parameters are 
increased, we take the free boundary condition in the two spatial directions. We found 
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FIG. 11: Spin correlation functions in / and // in phase diagram FigllOl Both Gn{r) and Gs{r) 
have no long-range order in / and //, however in // Gs{r) exhibits a short-range spiral order. 




FIG. 12: Spin correlation functions in /// and IV in phase diagram FigJlOl Gn{r) has long-range 
order both in /// and IV . In IV , Gs{r) exhibits a long-range spiral order, whereas it has usual 
long-range order in ///. 

that the correlators exhibit different behavior in the six "regions" I ~ VI shown in FiglTOl 
It is obvious that not only simple ferromagnetic (FM) correlation, antiferromagnetic(AF) 
correlations appear in these correlators. For example in the regions / and // (/// and 
IV), the correlation of the nematic order G„(r) exhibits the same behavior, but the spin 
correlator Gs{r) behaves differently in I and // (/// and IV), see FigUl] (Fig JT2i) . We 
have observed no phase boundary between the / and // (/// and IV) regions on which the 
internal energy E and specific heat C exhibit anomalous behavior. However from the result 
of the spin correlation function, Gs{r), we expect that the phases (u;^) ^ 0, {Or)) = and 
(u^) 7^ 0, {^{r)) 7^ are realized in the regions // and IV, respectively. On the other hand, 
there are no phases in the effective field theory that correspond to the phases V and VI of 
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the FM long-range order in the lattice model. 

Finally we turn on the C4-term in addition to the cs-term. Numerical study of the system 
was performed along the line C3 = C4 with ci = 6.0 and C2 = in the phase diagram. 
Observed C as a function of C3 = C4 is shown in FigJT^ There are two phase transition, 
one at 03 = 04^ 0.7 and the other at C3 = C4 ~ 6.0. In order to understand these phase 
transitions, we calculated the spin correlation functions in each phase, which are shown 
in FigJTSl From the result, it is obvious that the AF order disappears first and then at 
the second transition the spiral state appears. This state corresponds to the state with 
nonvanishing (^(r)) 7^ studied in the previous sections. 



VI. CONCLUSION 



In the present paper, we have studied the Sp(N) AF Heisenberg model in 2D that is 
expected to be realized in cold atom system in an optical lattice. We focus on the ground 
state structure of the system and derived the effective field theory for the system, which is a 
kind of extension of the CP^~^ nonlinear a-model. Then we studied the phase structure and 
critical behavior by using the 1/A^-expansion. We found that the spatial anisotropy induces 
a phase transition from the ordered state to paramagnetic state. As the explicit breaking 
of the SU(N) symmetry increases, the ground state exhibits a spiral order of the ad joint 
representation of Sp(N) group. This state also has a nematic order of vector representation 
of the Sp(N) group. 

We also introduced a lattice gauge-model counterpart of the effective field theory, and 




FIG. 13: Spin correlation functions in V and VI in phase diagram FigllOl Both Gn{r) and G's(r) 
exhibit the same behavior in V and VI. 
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FIG. 14: Specific heat on line C3 = C4 with ci = 6, C2 = 0. Result indicates two phase transitions 
on the hne. 




FIG. 15: Spin correlation functions for ci = 6.0 and C2 = 0. 

studied its phase structure by means of the MC simulations. We found a similar phase 
diagram with that of the field theory, but order of the phase transitions is different in two 
systems. 



In Ref. 13], finite-temperature phase diagram of Sp(N) model is studied by using a 
Ginzburg-Landau theory in terms of gauge-invariant spin fields, 0"* = z'^T'^^z, cf)"" = z'^T'^^z. 
There quantum phase transition of Sp(4) spin system in 3D stacked square lattice is also 
discussed. Phase transition from CP^ Neel ordered state to photon liquid state is predicted. 
In the (3+l)D counterpart of the lattice model (5.1), it is expected that deconfined photon 



phase exists 
previously 



br sufficiently large C2. U(l) gauge theory of CP^ field in (3+l)D was studied 
181], and it was found that phase transition from CP^ Neel ordered state to photon 
liquid state is of second order. 

It is interesting to study effects of hole doping to the Sp(N) AF magnets and investigate 
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how the long-range orders are broken and if a new phase with hole-pair condensation appears. 
This system is an extension of the t-J model for high-temperature superconductivity. This 
problem is under study and we hope that results will be reported in a future publication. 
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Appendix A: Derivation of the effective field theory 



In Sec. 2, we derived the effective field theory of the Sp(N) AF Heisenberg modeL In the 
present appendix, we give some details of the derivation. 

In Eg. (12. 141) . A{t) is expanded in powers of p and p. Vectors k and T are explicitly given 
as follows, 



v 



-U{woWi 

-U{W2WI 
-U{W3WI 

-U{WiW2 
-U{wiW3 



' ^ Jo,jizjJwo){zjJwi)U - ^ J'^j{woZj){zjWi)U ^ 
j j 

- Jo,j{ZjJ Wi){ZjJ Wi)U - J'^j{w2Zj){ZjWi)U 

j j 

- Y Jo,j{^jJ Wz){ZjJ Wi)U - Y J'o,jiwzZj){z.jWi)U 

i i / 

Y Jo,j{ZjJWo){ZjJWi)U - Y Jo,ji^j^o)i^iZj)U ^ 
j j 

Y Jo,j{ZjJW2){ZjJWi)U - Y J'o,MjW2){WlZj)U 

j j 
Yl ^3){zjJ wi)U - Y J'o,MjW^){wiZj)U 

i j I 



(A.l) 



(A.2) 



After integrating out p and p, we obtain Eq. (l2.17p . For J^^ ^ Jq^, the inverse of the matrix 
M can be approximated as 



M 



-1 

















AT . ~ 



iJo.i ' 4J, 



— T- 



-(4Jo,,)-i J 
^ 







(A.3) 
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and 



r ]^ 

odd k ^ 

- Jo,i 'ttt (^i^^fc) {ziJwi){wiWk) + (2;iJ'wfe)(ziJ^wi)(wfcWi)} 
k ^ i,j 

~ J'oi^' 7prY{{WkZi){ZiWi){wiWk) + {ZiWk){wiZi){WkWi)} 

k "^-^^ i 



+ {ZiJWk){ZiJwi){WkZj){ZjWi)} 



k i,J 



- Jo,i^Y{ZiJwi){ZiJwi) - j'^/Y{wiZi){ZiWi) 



(A.4) 



where Ylik' denotes the sum over k — 0,2,3, and Jq fa AJo,i + 4JpjAJ2 — Jz — 4:Jo,i. By 
using identity such as 



{ziZj) = {ziz){zzj) + a^{DizDjz), 
where a is the lattice spacing, and the completeness of {wfe}. 



(A.5) 



(A.6) 



fe=0 



we obtain 



odd 



1 T„2 



2Jd 



2d 



iZ 



^ ^{DrzDjZ - DjzDrz) + Ja^ ^ DjzDjZ - 2Jd 



'^2 



J'a 



^(zJDrz)izJDrz) + -i^Y.{-zJDfz){zJD,z) 



J'a^^{zJDjz){zJDj, 



Y^{{zJDjZ){zJDrz) - {zJDrZ){zJDjz)) 



2Jd 

+ ^ J2{izJDjz){z,JD^z) - {zJDjz)izjJD^z)} 



From FiglJl by substituting 



Dl = {-)yD, + Dy 

D2 = 2{-)yD^ 



the effective action Se defined by 



/3 



dTAJr) 



is obtained as 



/"^ r 1 - 

Se= dr —DrzDrZ + 2Ja^ V D^zD^z - a{-y{DrzD^z - D^zDrz) 
Jo A-J 



/i=x,y 



J' 
J' 



zJD,z){zJD,z) -2J'a^ ^ {zJD^z){zJD^i 
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J 
J' 



-)y{{zJD^z){zJDrz) - {zJDrz){zJD.,z)} 



+ — Y^{{zJD^z){z,JD^z) - {zJDjz){z^JD,z)] 
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By the rescaling I^M), 

S = — ; — / d^r D^,zDnZ — Vl + XadyiDrzDxZ — D^zDr^ 

\ / «^ l_ ,1— T- ^ 



- ^ V - ^Vl + Aa9j,(D,;2D,2 - D^zD^i 

Jo Jo 

^=T,x,y 
J' 

+ V -j-VTT\{{zJDjz)izjJD^z) - {zJDjz){zjJDrz)} 

i ^-^^ 

The last term in Eg. (I A. 11]) is rewritten as 

V -f Vl + \{{zJDjz){ZjJDrZ) - 
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Jo 



then, 



2A 



2a(l + A) 



where Sb is the Berry phase 

1 
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J d?r Yl D^zD^z-^^ Yl {zJD^z){zJD^ 
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(A.ll) 



(A.12) 



+ Sb, (A.13) 



(A.14) 
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